Why magnesium diboride is not described by anisotropic Ginzburg-Landau theory 
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It is well established that the superconductivity in the recently discovered superconducting com- 
pound MgB2 resides in the quasi-two-dimensional band (cr-band) and three-dimensional band (n- 
band). We demonstrate that, due to such band structure, the anisotropic Ginzburg-Landau theory 
practically does not have region of applicability, because gradient expansion in the c direction breaks 
down. In the case of dirty 7r-band we derive the simplest equations which describe properties of 
such superconductors near T c and explore some consequences of these equations. 



Ginzburg-Landau (GL) theory is the most powerful 
and widely used phenomenological theory of supercon- 
ductivity (see, e. g., Refs. GIB)- It describes practically 
all known superconductors in the vicinity of transition 
temperature. GL theory is fully microscopically justified 
and all its parameters can be derived from the micro- 
scopic BCS theory GL theory provides the basis for 
such elaborated fields as vortex physics Q and the theory 
of fluctuation phenomena . 

The recently discovered superconductor MgB2 gives 
an example of a superconductor which is not described 
by the anisotropic GL theory. This very unusual fea- 
ture is a consequence of a specific band structure of 
this compound. It is reliably established that super- 
conductivity in MgB 2 resides in two families of bands: 
strongly superconducting quasi-two-dimensional cr-bands 
and weakly superconducting three-dimensional 7r-bands 
(see, e.g., Ref. Both bands are characterized by 

their intrinsic coherence lengths, and the c axis coherence 
length in the cr-band is much smaller than c axis coher- 
ence length in the 7r-band. Typically, the strong band 
forces the order parameter in the weak band to change 
in the c direction at distances smaller than the intrin- 
sic c axis coherence length in this band. This means 
that almost in the whole temperature range the effec- 
tive coherence length in the c direction, £ Z (T), is smaller 
than the intrinsic coherence length in the 7r-band, £„. z . 
The crossover to the GL region takes place only when 
£ 2 (T) exceeds ^ T;Z , which occurs in the very close vicin- 
ity of T c . In this narrow region the 7r-band strongly in- 
creases the c axis coherence length. Beyond the narrow 
region, the variations of the order parameter in the c di- 
rections are not described by the anisotropic GL theory. 
Important consequencies of GL theory breakdown are 
the strong; temperature dependence of the H c2 anisotropy 
[3, la, 111 llfj and strong deviations of the H c2 angular de- 
pendence from the simple "effective mass" law h0. 111. 

Obviously, the breakdown of the anisotropic GL theory 
has numerous consequences and it would be desirable (i) 
to trace the reason of this breakdown and (ii) to derive 
the simplest model, which replaces the GL model near T c . 



This paper addresses these issues. For illustration, wc 
use the simplest microscopic model, multiband general- 
ization of the Usadel theory, describing a dirty two-band 
superconductor with weak interband scattering 0, Il2| . 
However, the main conclusions are very general and do 
not depend much on the intraband scattering strength. 
In the model we use the GL expansion for the cr-band 
and keep the microscopic description for the 7r-band, i.e., 
only " dirtyness" of the 7r-band is essential for a particular 
form of equations. 

We consider a dirty two-band superconductor with 
weak interband scattering. Such superconductor is de- 
scribed by Usadel equations for the impurity averaged 
normal and anomalous Green's functions, G a and F a , 



G~ 



■ 1, and the pair potentials A Q , 



V, 



uFct-Ysj [G a n)F a -F a V]Ga\ = A a G a (1) 



where a = 1,2 is the band index, j = x,y,z is the co- 
ordinate index, Dj = Vj — (2iri/&o)Aj, T> a ,j are diffu- 
sion constants, and oj = 2irT(s + 1/2) are Matsubara 
frequencies. Bearing in mind the application to MgB2, 
in our notations index 1 corresponds to cr-bands and in- 
dex 2 to 7r-bands, T>ij = T> a j and T> 2 j = T^nj- All 
bands are isotropic in the xy plane, T> a x — T> a ^ y and 
anisotropic in the xz plane with the anisotropy ratios 
7 Q = \/T>r l t/P» z- Self-consistency conditions can be 
written as 



a.x /'Da. 

m 



W 1 A 1 -W 12 A 2 = 2ttTJ2 ( F l " )+Ailni (2a.) 



t 



-IF 2 iA 1 + IF 2 A 2 = 27rT^ (f 2 -— )+A 2 lni (21,) 



t 



where t = T/T c and the matrix W a p is related to the 
matrix of coupling constants A. a p as 



-A_ + WA2+A 12 A2i k- + ^l\ 2 _+\ 12 k 21 

Wi = , W 2 - 



Det 

W 12 = Ai 2 /Det, W 21 = A 2 i/Det, 



Det 
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A_ = (An - A 22 )/2, Dct = A11A22 - Ai 2 A 2 i, W X W 2 
W\ 2 W 2 \. The supercurrent components are given by 



(3) 



a uj>0 



where N a are the partial densities of states 

We start with the derivation of the GL equations from 
the Usadel equations in the close vicinity of T c following 
a standard route. In the lowest approximation = 1 
and w A a /uj. When A Q are small and change slowly 
in space (the exact criterion will be derived below) one 
can keep only the leading nonlinear and gradient correc- 
tions 



U) 



A 3 T) 
2^ + ^j 2w* ' c 



(4) 



Substituting this expansion into the self-consistency con- 
ditions (J2J), we obtain coupled GL equations for two gap 
parameters |L3j 



W1A1 - W 12 A 2 



ei J D 2 J A 1 -bAl + rA 1 (5a) 
-W21A1 + W 2 A 2 = £ ZljD]A 2 - bA 3 2 + tA 2 (5b) 

with e aj = (n/8T)V ad , b = 7C(3)/(8tt 2 T 2 ), and r = 
ln(l/i) « (T c - T)/T c . 

Near T c the right hand sides of Eqs. (|H^,b) are small. 
This allows us to reduce Eqs. (0i,b) to a single GL equa- 
tion by looking for solution for A2 in the form 



A, 



W21 
W 2 



Ai + S 2 



(6) 



From Eqs. (|SJi,b) we obtain 

-w 12 s 2 = €i D2 j A i - 5A ? + rA i ( 7a ) 

W 2 5 2 =Y £ 1 D 2 A 2 -bAl + TA 2 . (7b) 



The second equation indicates that 8 2 is a small correc- 
tion, S 2 <C A2, and one can use A2 « (W21/W2) Ai in 
the right hand side of this equation. Excluding S 2 and 
introducing the band-averaged order parameter 

A 2 = W^l + WjA 2 _ W 12 Wj + W 2X Wl ^ 
W 2 + Wx ~ W 12 W 2 (W 2 + Wi) 1 

we finally obtain the anisotropic GL equation for A 

- ■ $ D j A + bA3 - rA = 0, (8) 
with the average coherence lengths 



W2 + Wi 



For the supercurrent, using relation W 2 \jW\ 2 = 
A21/A12 = Ni/N 2 , we derive 

jj w 4eN,P£ 2 Im [A*D.,-A] (9) 



with N = Ni + N 2 and 



a^ 2 (w 2 + Wx) 2 



{Ni + N 2 ) {N 2 W% + iVi W{) ' 



From Eqs. JSJ and J^J we derive the components of the 
London penetration depth 



A 



_ 2 32TT 2 eNP^r 



For the parameters of MgB 2 , Wi <C Wjjj £i. z "C £2,2, 
£13 ~ ^2,3;: the dominating effect of the 7r-band is the 
renormalization of the c axis lengths 



£l,z + ^12^2,2) 



A 



_, ^ 327r 2 erA^i 2 



^12^2,2) 



(10) 
(11) 



with 5i2 = Wi/W 2 < 1. The influence of the 7r-band 
on properties not related with the variations of the order 
parameter along the c axis are weak and can be treated 
perturbatively. 

We obtain now the applicability criterion for the 
GL expansion. The gradient expansion is justified if 
— £ 2 j V 2 A Q < A a for all a and i. Because a typical scale 
of the spatial variations is the temperature-dependent GL 
coherence length £j(T), this condition simply means 



&Cr)>£«j 



(12) 



The most restraining inequality is the one for a = 2 and 
i — z which gives 



(T c - T)/T c < gjZh + Si 



(13) 



Because £i j2 -C £2,2 and Si 2 <C 1, the applicability of the 
GL approach is limited to extremely narrow temperature 
range near T c , i.e., the situation is very different from 
conventional superconductors. For parameters of MgB2 
this condition implies (T c - T) /T c < 0.05. On the other 
hand, near T c the fluctuation effects become important. 
This means that the mean-field GL theory practically 
does not have a region of applicability. 

We derive now the simplest theory which replaces the 
GL theory in the conventional GL region (T c — T) /T c < 
1. As the gradient expansion actually breaks down only 
for the 7r-band, in the vicinity of T c we can proceed with 
the expansion @ for the c-band, a = 1. Substituting 
this expansion into the self-consistency conditions, we ob- 
tain Eq. (|5a|) . The 7r-band only weakly renormalizes the 
nonlinear term and we can use the linear approximation 
in this band 



^ 2 -£,%^ = a 2 . 



(14) 



The 7r-band order parameter can again be represented 
by Eq. © with 5 2 being a small correction. Finding 



3 



this correction from Eq. (|2b|l and substituting it into Eq. 
(|2a|l together with the GL expansion for Fi, we derive 
coupled equations for Ai and reduced 7r-band i^-function 
f s ,f s = (27rTW 2 /W 21 )F 2 (co s )M 

- (1 + 5i a ) rAi + 6A? - £ ^ ; />;A, 

-^E(/,-^J=0, d5a) 

(a + 1/2)/, - ^ ^. i J. = A x , (15b) 

and the expression for the supercurrent 

8 00 
./..- I'.Vi^ln, [AID i A 1 ]+-jiV 1 Si a £2 >J .£; Im [/^/J . 

These equations replace the GL equations in the case of 
dirty 7r-band. Note that the same equations are also valid 
in the case of clean a-band but with different definition 
of the coherence length £ ltj , ^ = 7C(3) /(4ttT) 2 . 

In the case of weak superconductivity in the 7r-band, 
S12 <C 1, and for £i )Z <C £2,*, one can neglect the in-plane 
gradients in the Eq. I|15b(l and obtain an even simpler set 
of equations which describe only the dominating strong 
effects, related to inhomogeneities of the gap parameter 
along the c axis, and neglect small renormalizations of 
the coefficients by the weak 7r-band 

-^E^-TTiVr ' (16a) 

(s + 1/2) f s - ^l z D 2 J s = Ai, (16b) 
./, « I'.V:-- , In, [A*D j A 1 ] 

+ s^rns^ E Im LfWJ > ( 16c ) 
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FIG. 1: Temperature dependence of the c axis coherence 
length, £ Z (T) computed from Eq. 1171 with parameters S12 = 
0.034 and „ = 300£? „. Marked GL region corresponds to 
condition £z(T) > £,r,z = £2,2- The inset shows dependence 
£z~ 2 (T) with the dashed line showing the linear GL asymp- 
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FIG. 2: Left panel: Dependencies j z (Pz) at different temper- 
atures for S r = 6. The curves are marked by the reduced 
temperatures t^2,z/Ci,z = 50, . . . , 300. In the unit of the ver- 
tical axis jdpi = 4eiVi£i ! z/& is the depairing current scaie for 
the a-band. Right panel: The temperature dependence of the 
depairing current for the same vaiue of S r - The dashed curve 
shows GL dependence. 



We explore now some consequences of these equations. 
To define an effective coherence length, we consider the 
response of the order parameter to the weak z-dependent 
variation of T c , r — > t(z) = t + 8t{z). In linear approxi- 
mation with respect to 8t(z) Eqs. i|16a|) and (|16b|) can be 
solved by Fourier transform yielding Ai = A^ + SAi(z) 

SA 1 (z) = J G(z-z')ST(z')dz', 

dk exp(ikz) 



G(z) 



2tt 2r + e iz k 2 + S 12 g[(2/^)e 2 z fc 2 ]' 



where g(u) = tp(l/2 + u) — ip(l/2) and ip(u) is the 
digamma-function. In contrast to the GL model, the 
decay of the perturbation 5A\(z) is not exponential. Us- 
ing the last equation, one can introduce the effective co- 
herence length which determines the scale of spatial 



variations of the order parameter in the z direction, 

ei,je z +s 129 [(2/n 2 )e 2 je z ]=r a?) 

The dependence £ Z (T) computed from this equation us- 
ing paran 
in Fig. 1. 

Consider the relation between the supercurrent j z and 
supermomentum p z = V z — (2ir/$o)A z , which deter- 
mines the c axis London length and depairing current. 
From Eqs. (|16af) and (|16b(l we obtain 



ipu t 

m 



ing parameters S12 = 0.034 [12| and £ 2 z — 300£ 2 z is shown 



jz(pz) = 4eiV 1 A 2 (p z )p z 



(18) 



x U?.+E 



V t^ (s + l/2+(2/^)(^ z p z ry 

A 2 (p 2 ) = (r- i\A - S 12 g [{2/^)e 2 , zPz }) lb. 
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In the linear regime j z w (AerNi/b) (£\ >z + Si2^, z ) Pz- 
This means that in the whole range (T c — T)/T c <C 1 
the z component of the London length is given by the 
GL formula (|1 1|) . In conventional superconductors the 
dependence j z (Pz) is nonmonotonic and its maximum 
gives the well-known GL result for depairing current, 
jd P = c<I> /(12\/37r 2 A 2 O oc r 3 / 2 for t -► 0. 2] In our 
case the situation is different. The amplitude of the or- 
der parameter is suppressed at p z ~ 1/£ Z (T). However, 
in the region £z(T) -C £2,2 the dependence j z (p z ) be- 
comes nonlinear at much smaller p z , p z ~ l/^2,z- The 
shape of this dependence is determined by the parameter 
S r = Si2^2. z /^i z - The dependencies j z {p z ) for S r = 6 
and different temperatures are plotted in the left panel 
of Fig. 2. For large values of S r the dependence j z (p z ) 
has two maxima within some temperature range, where 
first (second) maximum corresponds to the suppression 
of supercurrent in the 7r-(cr-)band. At low temperatures 
the depairing current jd P is given by the second maximum 
and is determined mainly by the cr-band. At certain tem- 
perature near T c global maximum switches to the first 
maximum (see Fig. 2). The temperature dependence of 
jd P has a kink at this temperature (see the right panel 
of Fig. 2). For S r > 6 the local maximum of j z (p z ) at 
Pz ~ 1/^2,2 exists even in the limit £*(T) <§C £2,2- 

As another example, we compute from Eqs. I jlfijl the 
in-plane upper critical field near T c , H c2ya {T). Exper- 
iment 0] shows strong upward curvature of H C 2, a (T), 
leading to the temperature dependent anisotropy factor. 
Microscopic calculations reproduce this feature, both in 
clean || and dirty cases, but require rather heavy nu- 
merical computations. Our model allows to trace origin 
of the upward curvature in a simple way. Selecting the 
gauge A z = Hx and introducing reduced variables h = 
H/H { £ with h£> = $o/(2<i,^ M ), x - Vhx/£ hx , 
r z = T>2 tZ /T>i^ z , we write the linear equation for deter- 
mination of the upper critical field, h = 



as 



c 00 



jfh^- fs )- wlAi+x2Ai=T h Au (19a) 



(s+l/2)f s + —r z hx 2 f s = A 1 . (19b) 



Excluding f s , we obtain the Schroedinger equation for 
Ai with non-parabolic potential 



2r z hx 2 



Ax = -Ar 
h 



(20) 



Only in the limit h -C vT + S\2r z /r z -C 1 this equa- 
tion reduces to the usual oscillator equation. In this 
limit, using expansion g(u) » (ir 2 /2)u, we reproduce 
the GL result, h C 2 = r/Vl + S\ 2 r z . The inequal- 
ity h C 2 -C \/l + Si2r z /r z reproduces criterion for 
the validity of the GL theory. In the opposite limit, 
VI + S\ 2 r z /r z <C h <C 1, one can use the asymptotics 



g(u) fn ln(u) + 7£ + 2 In 2 for u > 1 with 7^ w 0.577 
being the Euler constant, and obtain 

- V 2 Ai + fx 2 + ^ ln(a; 2 )") Aj - a(h)Ai 



t = a(h)h + Si 2 I In 



(1 „ 


~8hr z ~ 










. ^ 2 . 





This gives the following equation for the upper critical 
field 

h C 2 + S12 In [Cft, c2 r z ] = t 

with C ~ 1 (C sa (8/7r 2 )exp [(hi (re 2 )) + j E ] = 2/tt 2 
for h 3> S12 ). In this limit the 7r-band gives only small 
logarithmic correction to the upper critical field. As we 
can see, the upper critical field has a strong upward cur- 
vature in a narrow region neat T c : the slope dh c2 /dT 
changes from 1/V1 + Si2r z to 1 near r = S12 + l/r z , in 
agreement with microscopic calculations and experiment. 

In conclusion, we demonstrated that the properties of 
magnesium diboride are not described by the anisotropic 
GL theory. We derived a simple model which replaces 
this theory in the vicinity of T c , and explored some con- 
sequences of this model. 
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